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Abstract 



Our purpose in this paper is to study when a planar differential system 
polynomial in one variable linearizes in the sense that it has an inverse inte- 
grating factor which can be constructed by means of the solutions of linear 
differential equations. We give several families of differential systems which 
illustrate how the integrability of the system passes through the solutions of 
a linear differential equation. At the end of the work, we describe some fam- 
ilies of differential systems which are Darboux integrable and whose inverse 
integrating factor is constructed using the solutions of a second-order linear 
differential equation defining a family of orthogonal polynomials. 
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where P{x,y) and Q{x,y) belong to the ring of real polynomials in one variable 
and are analytic in the other variable, that is, they belong to the ring R(x)[|/] if we 
choose y as the variable in which they are polynomial. We will always assume that 
P{x, y) and Q{x, y) are coprime polynomials with respect to y. We denote by d the 
maximum of the degrees of P and Q as polynomials in y. 
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1 Introduction 



In this work we consider planar polynomial differential systems as: 




(1) 
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We define the orbital equation associated to system ([T]) as the ordinary differ- 
ential equation which is satisfied by the orbits of the system, that is, the orbital 
equation associated to system ([1]) is either 



dy ^ Q{x,y) 
dx P{x, y) ' 



or 



dx P{x,y) 
dy Q{x,y)' 



The aim of this work is to study when a system ([T]) linearizes. 

Definition 1 We say that system ([1]) linearizes, or that it is linearizahle, if it has 
an inverse integrating factor which can he constructed by means of the solutions of 
linear differential equations. 

We recall that the classical definition that a system ([T]) is linearizable is that 
there exists a change of variables which transforms the orbital equation associated 
to system ([1]) into a linear differential equation. The techniques used to find such a 
change usually come from the Lie group theory, see [3l [9] and the references therein. 
We do not treat this problem in this paper but the examples that we study show 
that there is a connection between both definitions of linearizability. 

This paper is related to the integrability problem which is defined as the prob- 
lem of finding a first integral for a planar differential system and determining the 
functional class it must belong to. We recall that a first integral H{x,y) of system 
([T]) is a function of class in some open set U of R^, non locally constant and 
which satisfies the following partial differential equation: 



An inverse integrating factor of system ([T]) is a function V{x,y) of class in some 
open set W of M^, non locally null and which satisfies the following partial differential 
equation: 



The function {dP/dx) + {dQ/dy) is called the divergence of system ([T]) and it 
is denoted by div throughout the rest of the paper. We note that the function 
1/V{x,y) is an integrating factor for system ([1]) in U, and that given an inverse 
integrating factor defined in U, a first integral in W — {V = 0} can be constructed 
by means of the following line integral: 



P{x,y)^{x,y) + Q{x,y)—{x,y) = 0. 
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where (xo,?/o) is any chosen base point in U with V^(xo,?/o) 7^ 0. We note that this 
function H{x,y) is only well-defined, in general, in a simply-connected subset of 
U — {V = 0}. In nonsimply-connected subsets of U — {V = 0}, H{x,y) can be a 
multivalued function but it continues to exhibit the dynamic behavior of the orbits 
in the set. 

The integrability of system ([T]) is given, in many occasions, by the existence of 
invariant curves. We say that a function / : W C ^ M is an invariant curve 
for a system ([1]) if it is not locally constant and satisfies 

df df 

with kf{x^ y) a polynomial in y of degree lower or equal than d — 1, where d is the 
degree of the system in and it is of class in the other variable. This function 
kf{x,y) is called the cof actor oi f{x,y). In case f{x,y) = defines a curve in the 
real plane, this definition implies that the vector field associated to system ([1]) is 
tangent to the curve / = 0. In case f{x, y) is a polynomial we say that /(x, = is 
an invariant algebraic curve for system ([T]). The construction of inverse integrating 
factors or analytic first integrals inside certain functional classes (polynomial, ra- 
tional, elementary or Liouvillian) is strongly related with the existence of invariant 
algebraic curves, see for instance the work ^ and specially the references therein, 
and it belongs to the context of the Darboux theory of integrability. 

When considering the integrability problem we are also addressed to study how 
the existence of a first integral in a certain functional class implies the existence of 
an inverse integrating factor inside a certain given class of functions. In the partic- 
ular case that system ([1]) is polynomial, we have that the existence of an elementary 
first integral implies the existence of an inverse integrating factor which is a rational 
function up to a rational power. Moreover, when system ([T]) is a polynomial, we 
have that the existence of a Liouvillian first integral implies the existence of an 
inverse integrating factor of Darboux type, see [6], |8] and the references therein for 
the proof of these results. A function of the form exp [h/g) fi^f2^ . . . fg% where 
h{x, y) and g{x, y) are polynomials, each y) = is an invariant algebraic curve 
for system ([I]), Aj G C not all of them null, for i = 1, 2, . . . , s and s G N, is called a 
Darboux function. These results suggest that the functional class of an inverse in- 
tegrating factor is usually easier than the functional class of a first integral. This is 
the reason why we look for an inverse integrating factor to study the integrability of 
system ([1]). Moreover, the inverse integrating factor is shown to be defined in phase 
portraits in which the dynamics avoid the existence of a first integral. However, 
there are also systems whose dynamics avoid the existence of an inverse integrating 
factor, see [5j. 

In the work [8], systems of the form ([1]) whose integrability is given by the so- 
lutions of linear differential equations are described. We obtained a result which 
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allows to find an explicit expression for a first integral of a certain type. By means 
of a rational change of variable, we made correspond the homogenous second or- 
der linear differential equation: A2{x)w" (x) + Ai{x)w'{x) + Aq{x)w{x) = 0, whose 
coefficients are polynomials, to a planar polynomial differential system. We prove 
that this system has an invariant curve for each arbitrary non null solution w{x) of 
the second-order ordinary differential equation, which, in case w{x) is a polynomial, 
gives rise to an invariant algebraic curve. In addition, we give an explicit expression 
of a first integral for the system constructed from two independent solutions of the 
second order ordinary differential equation. This first integral is not, in general, a 
Liouvillian function. The inverse integrating factor of the system ([1]) which is associ- 
ated to the aforementioned second order linear differential equation, takes the form 
V{x,y) = q{x){w'{x) — g{x,y)w{x)Y where q{x) = A2(x)exp|J Ai{x) / A2{x) dx^ , 
g{x,y) is a fixed rational function and w{x) is a non-null solution of the second 
order linear differential equation. 

Moreover, in the work |8] we also consider first order linear differential equations: 
Ai{x)w'{x) + Ao{x)w{x) = with polynomial coefficients and analogous results are 
obtained. The inverse integrating factor of the system ([T]) which is associated to this 
first order linear differential equation, takes the form V{x, y) = Ai{x)g{x, y){w{x) — 
a{x,y)) where a{x,y) is a function defined in terms of Ai{x) and Ao{x), g{x,y) is 
fixed rational functions and w{x) is a non-null solution of the first order linear 
differential equation. 

Hence, in the work [5], we give families of systems which, by construction, lin- 
earize, because its corresponding inverse integrating factor is obtained in terms 
of the solutions of a linear differential equation. Moreover, the given families are 
very general since they come from any rational change of variables. The present 
work arises as a reciprocal of the work [8], since we look for systems which can be 
linearized, in the sense of Definition [H 

The goal of this work is to demonstrate an algorithm to detect when a system 
is integrable (either inside the Liouvillian class or not) by means of a linearization 
process. That is, we target to find systems whose integrability passes through the 
solutions of a linear differential equation. We proceed by giving and explaining 
several examples which illustrate this process. 

The examples that we study suggest that the integrability by linearization of a 
polynomial system ([1]) reduces to solve linear differential equations of order 2 or it 
falls into the Darboux theory of integrability. 

The studied examples also make arise the following questions: when a system 
is linearizable (in the sense of Definition [1]) with a linear differential equation of 
order 1, does it always exist a rational change of variables which transforms the 
system to an orbital equation which is linear? In the same way, we can also ask 
whether when a system is linearizable (in the sense of Definition [1]) with a linear 
differential equation of order 2, does it always exist a rational change of variables 
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which transforms the system to an orbital equation which is of Riccati type? 

For the famihes of systems studied in [H], the answer to the previous two ques- 
tions is affirmative. 

The question of hnearizabihty has attracted many authors since the transfor- 
mation of an ordinary differential equation or a partial differential equation of any 
order by means of several differential-algebraic manipulations to a linear differential 
equation, gives in general the solution of the first, nonlinear problem. Moreover, 
ordinary differential equations which linearize come naturally in some physical ap- 
plications, see [7] and the references therein. In the work [7] , the question of which 
ordinary differential equations (of any order) linearize upon differentiation is ad- 
dressed and some sufficient conditions on the form of the equation are given. How- 
ever, these sufficient conditions are very restrictive over the equations and only very 
special particular equations can satisfy them. We only consider ordinary differen- 
tial equations of first order, that is systems of the form ([1]) and we study several 
differential-algebraic manipulations so as to get a linear equation which character- 
izes its integrability. 

We use two different methods to exhibit that a system is linearizable: equivalence 
and compatibility. Both methods start in the same way. We consider a system ([1]) 
and we think of it as polynomial in one variable, for instance y. Then, we take a 
polynomial in the variable y of a certain fixed degree and with arbitrary coefficients, 
which are functions of the variable x, YlfLo hi{x)y\ and we impose it to be an 
inverse integrating factor of the corresponding system ([1]). This condition gives rise 
to a system of linear differential equations on the coefficients hi{x). In general, 
this system of linear differential equations is over determined. Several conditions on 
system ([T]) can make this system compatible and the way to choose these conditions 
is what distinguishes between both methods. 

The compatibility method is the most gross: we consider the system of linear 
differential equations with variables hi{x) and we uncouple the variables by means of 
differentiation and resultants. We end up with an algebraic-differential condition on 
system ([T]). Although this method gives all the possible choices of system ([T]) to have 
an inverse integrating factor of the prescribed form, it is usually too overwhelming 
to be carried out. 

On the other hand, we can consider the equivalence method. This method is wiser 
and consists on uncoupling the system of linear differential equations with variables 
hi{x), trying to avoid differentiation and resultants and only using substitution at 
each step. We end up with a number of linear differential equations of certain order 
and of only one variable, say ho{x), and we make these equations equivalent, that 
is, we impose them to be the same equation and/or to be identically null for some 
of them. This way gives certain particular conditions on the system which are. 
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usually, easy to satisfy. 

We note that the conditions given by the equivalence method are also contained 
in the conditions given by the compatibility method but their determination is much 
easier when equivalence is involved. 

The method of equivalence gives rise to one linear ordinary differential equation 
of a certain order ^ for one variable, which can be any hi{x) in the expression 
/ij(a;)?/\ Each solution of this ordinary differential equation gives, by back 
substitution, an inverse integrating factor for system ([1]). We would like to know 
the possible values of I, that is, we ask whether we can linearize systems ([1]) by 
means of a linear ordinary differential equation of any order £, with £ > 0. We aim 
to know the values of the order C. corresponding to a linear differential equation such 
that each of its solutions cannot be expressed as a polynomial on the solutions of an 
equation of lower order. The following result shows that the order of such a linear 
differential equation is at most 2. 

Theorem 2 We assume that system ([T]) has an inverse integrating factor of the 



where M is a nonnegative integer number and hi{x) are analytic functions in x, 
i = 0, 1, . . . , M . We assume that the functions hi{x), for i = 1, 2, . . . , M , are poly- 
nomials in ho{x) and its derivatives, and that ho{x) satisfies a linear differential 
equation of order i, with £ a nonnegative integer, whose solutions cannot be alge- 
braically expressed in terms of the solutions of an equation of lower order. Then, 
£<2. 

Proof. We know that a linear ordinary differential equation of order i has a fun- 
damental set of solutions with cardinal i. That is, there are i linearly indepen- 
dent solutions to the equation. Assume that i > 3 and let Vi{x,y), V2{x,y) and 
V3{x,y) be three inverse integrating factors each one constructed by using one of 
these linearly independent solutions through the expression ([2]). The quotients of 
two of them give first integrals of system ([T]): Hi{x,y) = Vi{x,y)/V3{x,y) and 
H2{x,y) = V2{x,y)/Vi{x,y). These two first integrals need to be functionally de- 
pendent since any first integral of a planar differential system like ([T]) is a function 
of another one. We are going to show that, in fact. Hi and H2 are algebraically de- 
pendent, that is, there exists a polynomial with real coefficients P{zi, Z2) such that 
P{Hi,H2) = 0. We consider the level curves of each Hf Vi{x,y) — CjV3(x, y) = 0, 
with i = 1,2, which are two polynomials in y because each Vi{x, y), i = 1,2, 3, is a 
polynomial in y. Let us take the resultant of the polynomials Vi{x,y) — CiV3{x,y) 
and V2{x, y) — C2Vz{x, y) with respect to y and we denote it by R{ci, 02, x). We re- 
mark that this resultant is a polynomial in ci and C2 and we are going to show that 



form: 



M 




(2) 



i=0 



6 



it factorizes as R{ci,C2,x) = P{ci,C2)Ro{ci,C2,x) where P(ci,C2) and Ro{ci,C2,x) 
are polynomials in ci and C2- This factorization of R{ci, C2, x) is deduced by the fact 
that each y-root of Vi{x, y) — ciV3(a;, y) = 0, of a fixed ci, needs to correspond to a 
value of C2 such that the whole y-root is contained in V2(x, y) —C2V^{x, y) = 0. That 
is, fixed a ci and a y-root of Vi{x, y) — ciV^^x, y) = 0, there exists a value of C2 such 
that this y-root is completely contained in V2{x, y) — C2V3(x, y) = 0. Therefore, we 
have encountered a polynomial P(ci, C2) which relates the two first integrals in the 
desired way. 

Let us call Si{x) the solution of the linear ordinary differential equation which 
gives the inverse integrating factor Vi{x,y), i = 1,2,3. Since P{Hi,H2) = 0, we 
deduce that there exists a homogeneous polynomial with real coefficients such that 
p{Si, S2, S3) = 0. To encounter this polynomial p, we put y = in the expression of 
P{Hi, H2) and we take a common denominator. The existence of this polynomial 
p implies that 6*3 can be algebraically expressed in terms of 5*1 and 82- We remark 
that any two given functions and 5*2 satisfy the linear homogeneous ordinary 
differential equation of second order: 



det 



w"{x) w'{x) w{x) 
S'lix) S[ix) S,{x) 
S'^ix) S',{x) S2{x) 



0. 



Thus, the function 5*3 is algebraically expressed in terms of the solutions of an 
equation of second order. 

We conclude that any inverse integrating factor which shows the linearizability 
of a system ([T]) through a linear differential equation of order i with i > 3, can be 
expressed in a way that the linearizability of the system is given through a linear 
differential equation of, at most, second order. | 

In order to make more precise the given notion of linearization of a system ([T]), 
we include several examples of this phenomenon and we use the described methods 
of equivalence or compatibility. 



2 Examples 

2.1 Automatically linearizable systems 

In this section we describe examples of systems for which we find an inverse integrat- 
ing factor constructed by means of the solutions of a linear differential equation, 
that is, we describe examples of linearizability. We do not need to impose any 
condition on the system to ensure its linearizability, that is, the system of linear 
differential equations on the functions hi{x) for which V{x,y) = J^fLo hi{x)y^ is an 
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inverse integrating factor is not overdetermined. 

Example 1. Let us consider systems of the form ([1]) with: 

m 

P{x,y) = -y, Q{x,y) = Y,92^{x)y^\ (3) 

1=0 

where g2i{x) are analytic functions and m is an integer number with m > 0. In this 
section we only consider the case in which m = 2. 

We look for an inverse integrating factor V{x, y) which is a polynomial in y of 
the same degree as ([3]) and of the form: 

m 

V{x,y) = J2hMy'\ 

1=0 

where h2i{x) are suitable functions which will satisfy linear differential equations. 
Our goal is to impose such a function V{x,y) as an inverse integrating factor for 
system ([3]) and to deduce the relations on the functions g2i{x) to accomplish it. 

The case m = is easily integrable since the corresponding orbital equation has 
separated variables. Let us explicit the computations made when m = 2. We have 
the system: 

X = -y, y = go{.x) + g2{x)y^ + gi{x)y^, (4) 

and we look for an inverse integrating factor of the form V{x, y) = ho{x) + h2{x)y'^ + 
h4{x)y'^. We impose the following relation to be satisfied: 

dV dV 
ox ay 

where div = 2g2{x)y + 4:g4{x)y^. The previous partial differential equation can be 
arranged in powers of y and we get that three relations among the functions h2i{x) 
need to be satisfied, corresponding to the powers y^, and y. We remark that we 
have three relations and three functions to satisfy them. This is an exceptional case 
since the number of equations and the number of variables coincide. These relations 
read for: 

h'^{x) - 2g2{x)hi{x) + 2gi{x)h2{x) = 0, 

h'^ix) ~ 4:go{x)hi{x) + 4:gi{x)ho{x) = 0, (5) 

h'^i^x) - 2go{x)h2ix) + 2g2{x)ho{x) = 0. 

We can deduce the values of the functions h2{x) and /io(x), for instance, from 
the first two equations ([5]). The third equation in ([5]) gives a third order hnear 
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differential equation for /i4(x) which is: 



glh'H - 35(45-4/14 + {-^glgl + 165ofi'4 - 45452 + 452545I + - 9i9'i)K + 



where we avoid to write the dependence on the variable x to simplify the notation. 

We have obtained a third order linear differential equation whose solutions give 
rise to an inverse integrating factor for system ([3]) with m = 2. We note that 
we did not need to impose any restriction on the functions g2i{x), i = 0, 1,2, so 
as to linearize the system. System ([3]) with m = 2, that is system (jl]), is always 
linearizable and we have been able to deduce this fact using our method. 

Although we have obtained a third order linear differential equation for h^lx), 
we are going to see that this function is related to a second order linear ordinary 
differential equation, as the statement of Theorem [2] establishes. Let us impose an 
inverse integrating factor for system (jl]) of the form V{x,y) = (ho^x) + h2{x)y'^)^, 
where ho{x) and h2{x) are suitable functions. Repeating the same computations as 
before, we get that the relation corresponding to the power in the equation of 
inverse integrating factor gives that ho{x) = {g2{x)h2{x) — h^i^x)) / {2gi{x)) and the 
other two relations (corresponding to y^ and y^) are equal and give the following 
second order ordinary differential equation for h2{x): 



The fact that the two relations, corresponding to y^ and y^, are equal is not expected 
and this equality confirms that system (jl]) is linearizable. 

Actually, straightforward computations show that if we denote by A{x) and B{x) 
two independent solutions of the second order equation (j7]), then a fundamental set 
of solutions of (jH]) is: A'^{x), A{x)B{x) and B'^{x). Hence, the third order linear 
ordinary differential equation (jH]) is, in fact, reducible to a second order ordinary 
differential equation. This is an example of the result stated in Theorem [2l We note 
that, in any case, we have linearized system (jlj). The reduction of order passes, in 
this case, through the change h^lx) = {h2{x)Y which transforms equation (jOj) in a 
nonlinear differential equation which is compatible with (j7|). 

We need to make precise that when we have linearized system (j3]) with m = 2, 
that is system (jlj), we have obtained a third order linear differential equation (j6|), but 
this equation was not necessary since an inverse integrating factor can be obtained 
through a second order linear ordinary differential equation, as Theorem [2] states. 

We consider the algebraic change of variables y ^ z with z = ho{x) + h2{x)y'^, 
where hi, i = 0,2 are the functions which define the encountered inverse integrating 




(6) 



+ (252 + ^92)9^94 - ^929'i + g2gigl]hA = 



54 h'^ - 54 h'^ + (5254 - 929i + 4506-4 - 9l9i) h2 = 0. 



(7) 
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factor V{x,y) = {Iiq^x) + h2{x)y'^y. This change of variables apphed to system (jlj) 
gives the following orbital equation: 

dx h2{x) 

which is of Riccati type. 

An algebraic change applied to system (jlj) which gives an orbital equation of 
Riccati type is not unique. We note that the algebraic change of variables y ^ u 
with u = y"^ applied to system (jl]) gives the following orbital equation: 

— = -2{go{x) + g2ix)u + gi{x)u'^), (8) 

which is also of Riccati type. The linearizable systems studied in the examples show 
the existence of a rational change of variables which transforms the system to an 
orbital equation of Riccati type, as a counterpart. We do not look for changes of 
variables but for linearizability. 

It is well-known the equivalence between ordinary differential equations of Ric- 
cati type and second order linear differential equations. For instance, by the change 
u w with dw/dx = 2g4{x) u{x) w we have that the Riccati equation ([8]) is equiv- 
alent to the following second order linear differential equation for w{x): 

g4,{x)w"{x) + {2g2{x)g4,{x) - g'^{x))w'{x) + 4:go{x)gl{x) w{x) = 0. (9) 

This family of systems also appears as a particular case of the systems described 
in [5] by means of the algebraic change of variables y ^ u with u = y"^, which leads it 
to the orbital equation of Riccati type ([H]). Following the ideas described in [S], if we 
consider equation ([7]) and we perform the change h2{x) = exp [J^ g2{s)ds] w{x), 
we obtain the linear differential equation (Q. We would like to remark that the 
linearizability process does not pass through the second order linear differential 
equation Q, which is equivalent to the orbital equation (jH]). The linearizability of 
system is concerned with the second order linear differential equation ([7]) whose 
solutions define an inverse integrating factor for the system. However, to impose 
that system (jl]) has an inverse integrating factor constructed with the solutions of 
a linear differential equation seems to imply that the orbital equation associated to 
the system is equivalent to a linear differential equation by means of an algebraic 
change of variables. All the examples that we present in this work confirm this 
implication although the involved linear differential equations are different and come 
from different sources. 

2.2 Linearizability by equivalence 

In this section we describe several examples of systems of the form ([1]) which, under 
certain restrictions, are linearizable. The way we determine these restrictions is by 
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using the equivalence method. 



Example 2. We consider system ([3]) with m > 3 and we explain a process 
which encounters a linearizable subfamily. Let us describe the computations for 
system ([3]) with m = 3. In fact, for m = 1 and for any m > 3 the discussion is 
analogous. The only cases which are different are m = (separated variables) and 
m = 2 (second order linear differential equation) which have already been treated 
in Example 1. We consider the system: 



X = -y, y = 9g{x) + g2{x)y'^ + gi{x)y^ + ge{x)y 



and we take a function of the form V"(x, y) = ho{x) + h2{x)y'^ + h4^{x)y^ + hQ{x)y^. 
By imposing it to be an inverse integrating factor and equating the same powers 
of y, we get 5 relations which need to be satisfied and which correspond to the 
coefficients of y, y^, y^, y'^ and y'^ . From the coefficient of y'^ we compute h/i{x) in 
terms of g2i{x) and he{x). In the same way, from the coefficient of y'^ we compute 
h2{x) and from the coefficient of y^, we compute ho^x) in terms of g2i{x) and Hq^x). 
We are left with two linear differential equations of second order for hQ{x) which 
read for: 



4^6 v3^i ^gl 



1 ( 9l\ glg'e f 92 



^96 \96j ^9i Vfl'6 



h = 0, 



94 K ^ l^^y _ _94_ ^ 9294 _ ^ ] ^/ ^ 



12gQ \3 ygaj Ug^ Gg^ 2 



, ( 949'i 9496 ,92f94\. 9o96 ;\ , _ n 

We apply the equivalence method to these two linear differential equations, that is, 
we impose the values of g2i{x) to make them the same equation. Astonishingly, we 
only need to impose the condition: 



x'^ 




, ^ 92{x)g4{x) 2 5(41 

9q{^) = o / A 

3gQ{x) 

so as to get only one second order linear differential equation for he^x): 

3gX + i^gl-l2g2g6-H)K + {^g^g', - ^ + I2g2g'e - I2geg'2^ h = 0. (10) 

Hence, we have a family of systems of the form ([3]) which linearize. At this point 
we have met our target since we have encountered a family of systems whose inte- 
grability passes through the solution of a second order linear differential equation. 
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We are going to describe another unexpected phenomenon which occurs in this 
family. The second order hnear equation for hg can be reduced to a hnear equation 
of order 1. We remark that Theorem [2] does not apply because equation ffTOj) is of 
order 2. We consider system ([3]) with m = 3 and the described value of the function 
go{x). We impose a function of the form V{x,y) = {ho{x) + h2{x)y^Y to be an 
inverse integrating factor. We get that Hq = g4:h2/{3gQ) and we only obtain one 
linear homogeneous differential equation of order 1 for the function h2{x) which is 
9(76^2 + 4(5'! — 3g2g&)h2 = 0. Therefore, we have that system ([3]) with m = 3 and 
the described value of go{x) linearizes. We remark that the final equation for h2{x) 
in the case m = 3 is of order 1 whereas the final equation for h2{x) in the case 
m = 2 is of order 2. 

Moreover, this family of systems also appears as a particular case of the work 
fSj. We consider the algebraic change of variables y ^ z with z = ho{x) + h2{x)y'^, 
where hi, i = 0,2 are the functions which define the encountered inverse integrating 
factor V{x, y) = {ho{x) + h2{x)y'^Y . This change of variables applied the considered 
family of systems gives the following orbital equation: 

dz 2 / g^jxf 9geix)z^\ 

dx 9 V 9e{x) h2{xf J 

which is a differential equation of Bernoulli type. The change z ^ u with z = 1/ a/u 
transforms equation ( ITTi) to a linear differential equation. We note that by this 
method we impose a system to be linearized, that is, the other way round of what 
we obtained in [8], where we started by the linear differential equation and we 
deduced the corresponding system. As far as the examples described in this work 
indicate, we have observed that when a linearizable system ([1]) (in the sense of 
Definition [T]) admits an inverse integrating factor of the form V{x,y) = c{x,y)"', 
where n in a positive integer and c{x, y) is a polynomial in y, the algebraic change 
y ^ z where z = c{x, y) applied to the system gives an orbital equation which is of 
Riccati or Bernoulli type. 

The same computations can be done for the systems ([3]) with m = 1 or m > 
3, and giving certain values to m — 2 of the functions g2i{x), that is we have 3 
arbitrary functions g2i{x) in the expression of Q{x,y), an inverse integrating factor 
of the form V{x,y) = {ho{x) + h2{x)y'^)"^ is exhibited, where ho{x) is expressed in 
terms of the g2i{x) and h2{x), and the function h2{x) is the solution of a linear 
homogeneous ordinary differential equation of order 1. In the particular case where 
the linearizable system with m > 3 is a polynomial differential system, we have 
that the encountered inverse integrating factor is of Darboux type. We remark that 
in the method of linearization we are not seeking for invariant algebraic curves but 
we have obtained them by construction. 

When we have studied systems of the form ([3]) with m = 2 we have exhibited a 
system in which the linearization is automatically met. When we have described the 
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systems of the form ([3D with m = 3, we have given a family of hnearizable systems 
in which the hnearization is met by the equivalence method. We also remark that in 
case system (|3]) is polynomial, we can encounter non Liouvillian inverse integrating 
factors when m = 2, that is, when the system is of degree 4 in y. When m > 3 in a 
polynomial system ([3]), the linearization process gives a Darboux inverse integrating 
factor and, thus, a Liouvillian first integral. This fact implies that the encountered 
second order linear differential equation flTUl) only has elementary solutions. 

Example 3. The following family of systems is quite more general and we are 
going to use equivalence so as to encounter a subfamily of hnearizable systems. Let 
us consider the system: 

X = co(x) + C2(a;)y^ 

y = do{x) + di{x)y + d2{x)y^ + d3{x)y^ + d4{x)y'^ + d5{x)y^ + d(i{x)y^ , 

where Ci{x) and di{x) are analytic functions with C2(x) ■ d^^x) ^ 0. 
We start by imposing an inverse integrating factor of the form 



i=0 



which is a polynomial in y of the same degree as system f[T^ . We substitute this 
expression of V in the partial differential equation that must be satisfied to be an 
inverse integrating factor and we obtain that a certain polynomial in y of degree 10 
must be identically zero. From the coefficients in y of this equation of degrees from 
10 to 5 we deduce the values of hi{x) with i = 0,1, ... ,5. The rest of the coefficients 
give five linear differential equations for h^lx). Two of them are of order 3 and the 
rest are of order 2. We impose the three equations of order 2 to be identically null, 
which give the following conditions on system (|T2l) : 



6co(x)4(x) 9co(x)24(x) 
d^{x) = 0, d4(x) = — , d2{x) - 



C2{x) ' C2{xy 

3co(x)rf3(x) / co(x) 

di[X) = — — — C2{X) 



C2{x) \C2{x) 



(13) 



Under these conditions we get that the three equations of order 2 are identically null 
and, surprisingly, the two equations of third order are equal. The fact that under 
this small number of restrictions, the involved linear differential equations become 
equivalent is unexpected and confirms the hidden structure of the linearizability 
process. This subfamily of systems (1121) is hnearizable. 

We remark that this is not the only way to proceed so as to get linearizability. 
We could also have imposed the two equations of third order to be identically null 
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and then use equivalence with the second order equations. We have only presented 
one of the possible cases that can be encountered. 

We know that the considered third order linear differential equation can be 
reduced because its solutions cannot be functionally independent. Let us consider 
an inverse integrating factor for system f|T2|) of the form V{x,y) = {h3{x)y^ + 
h2{x)y'^ + hi{x)y + ho{x)Y. Analogous computations give that, under the same 
conditions, the determination of V{x, y) comes from the solution of the following 
second order linear differential equation in h^{x): 

We note that the change of variables y ^ u with 

C2{xYu + 3(i6(a;)?/(3co(x) + C2{x)y'^) = 

transforms the orbital equation associated to system ( |T2|) . with the values imposed 
in (fT3|) . to a Riccati equation. We have that system (fT2|) . with the values described 
in (fT3l) . is a particular case of the results given in [8]. As before, when we apply the 
linearization process we do not look for systems which come from Riccati equations 
via a change of variables, but we obtain such systems as a counterpart. 

2.3 Linearizability by compatibility 

In this section we describe several examples of systems of the form ([1]) which, under 
certain restrictions, are linearizable, in the sense of Definition [H We encounter these 
restrictions by applying the compatibility method. 

Example 4. Let us consider the following system which appears in the work 
[Ij. The system 

X = y{-l + 2p^{x^ - y^)), y = x + px"^ + py'^ + 4p^xy^ (14) 

where p G M, has the inverse integrating factor V{x,y) = (x^ + y^)^(l + 2px + 
p^(x^ + y"^))- We are going to re-encounter this inverse integrating factor using 
compatibility. This is the first example in which the compatibility method appears. 

We consider system f ll4p and we impose an inverse integrating factor of the form 
V{x, y) = Hq^x) + h2{x)y'^ + h/i{x)y'^ + hQ{x)y^. We equate to zero the coefficients of 
the powers of y in the relation that makes V{x, y) an inverse integrating factor. We 
have five equations corresponding to the coefficients of y^ for j = 1, 3, 5, 7, 9 which 
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read for: 

2x{l + px)h2{x) = (1 - 2p^x^)/io(a;) + 2p(l + &px)ho{x), 

4a;(l + px)hi{x) - 2p^h'^^{x) = (1 - 2pV)/i2(a;) + Ap^xh2{x), 
6x(l + px)hf,{x) - 2p^h'^{x) = (1 - 2p^x^)h'^{x) - 2p(l + 2px)hA{x), 
2p^h'^{x) = -{l-2p^x^)h'e{x)+4p{l + 3px)he{x), 
2p^he{x) = 0. 

Wc remark that wc have five hnear differential equation for the four functions hQ{x), 
h2{x), h^i^x) and hQ{x). This system of hnear differential equations is shown to be 
overdetermined. We can proceed in several ways but we do not loss any generality in 
following one of them. We can, for instance, take the equations from the last to the 
first one in the order they have been written. We solve them, leaving an arbitrary 
constant at each step. We end up with several algebraic relations for these constants 
which mark their value and from which we get the previously described polynomial 
inverse integrating factor. The compatibility of these relations is obtained by an 
adequate choice of the constants of integration. Another way to study this system 
of linear differential equations is to take them in the order they have been written. 
From the first one we equate h2{x) and we substitute in the rest of equations. From 
the second equation, we equate h,i{x) and we substitute in the rest of equations 
and from the third equation we get Hq^x). We end up with two linear differential 
equations of fourth order for ho{x) and we make them compatible. The compatibility 
process goes as follows: we consider the two linear differential equations of fourth 
order for ho{x) and we make a linear combination of them so as to get a linear 
differential equation of third order. We derive this third order linear differential 
equation and we combine it with one of the previously considered linear differential 
equations of fourth order for ho{x). We have two third order linear differential 
equations for ho{x) and we combine them so as to get a linear differential equation 
of second order. We derive it and we combine with one of the previously considered 
equations of third order. We have at this step two linear differential equations of 
second order for ho{x), which we combine so as to get a first order linear differential 
equation for ho{x). We derive it and we obtain a second order linear differential 
equation for ho{x) which, combined with one of the previous gives rise to a first 
order linear differential equation. The two first order linear differential equations 
for ho{x) turn out to be the same. If this was not the case, we would combine them 
and we would obtain a compatibility condition on the coefficients of the system. In 
our case, we solve this first order linear differential equation for hQ{x). This value 
of ho{x) is the value which makes compatible the two fourth order linear differential 
equations from which we started the process. The only possible common solution of 
these two fourth order linear differential equations is ho{x) = x'^{l + px^ (modulus 
a multiplicative constant). We observe that this ho{x) univocally determines the 
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previously described inverse integrating factor. 

This example suggests that the integrability by linearization of a polynomial 
system reduces to solve linear differential equations of order 2 or it falls into the 
Darboux theory of integrability. We remark that any Darboux inverse integrating 
factor which is a polynomial in the variable y is encountered by our linearization 
process: either by the equivalence method, either by the compatibility method. 

We note that when applying the linearization process to system f[T^ we obtain 
a Darboux inverse integrating factor and, thus, invariant algebraic curves of the 
systems as a counterpart. 

Example 5. In this example we are addressed the question whether the family 
f lT^ can be embedded in a linearizable family and it can be seen as a particular 
case of linearization with equivalence. 

We can think of system f|T^ as a particular case of the following family of 
systems: 

X = y (giix) + g2ix)y^) , y = g^ix) + gA{x)y'^, (15) 

where the gi{x) are arbitrary functions. In case g2{x) = 0, we have that system (1151) 
coincides with system ([2]) with m = 2 after a time-rescaling, which has already been 
studied in the first example. Therefore, we can assume, without loss of generality, 
that g2{x) = 1. We note that the family f|T5|) is reversible by the change {x,y,t) \—>- 
(x, —y, —t) and, thus, any ansatz of an inverse integrating factor must be even in y. 
Let us consider an inverse integrating factor of the form: 

V{x,y) = ho{x) + h2{x)y'^ + hi{x)y^ + hQ{x)y^. 

The imposition for it to be an inverse integrating factor gives a polynomial in y. We 
observe that this relation is an odd polynomial in y of degree 9 and which is odd. 
The functions hi{x) must vanish each one of the coefficients of this polynomial in y. 
From the coefficient of y^ of this polynomial we deduce that h^lx) = k^, with a 
constant value which we assume to be nonzero. From the coefficient of y, we deduce 
the value of h2{x) and from the coefficient of y^, we deduce the value of /i4(x). The 
coefficients of y^ and y"^ give rise to two linear differential equations of order 3 for 
ho{x). 

When applying equivalence to these two equations, that is, imposing them to be 
the same equation, we deduce the following conditions: 

1 , g' (x) 

gsix) = -{2gi{x)g4{x) - gi{x)g[{x)) , g^ix) = - 

In this case, the orbital equation associated to system ( ITSl) is of separated variables. 
We observe that the linearization process for this example leads to an ordinary 
differential equation with separated variables. 
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Another way to study the possible hnearizabihty of system fllSI) is to impose 
that the two hnear equations of third order have a nonzero common solution, that 
is, to impose compatibility. To perform the computations of compatibility for two 
equations of third order carries long calculations and many cases. If we apply com- 
patibility to the two linear differential equations of order 3 for ho{x), we obtain sev- 
eral conditions. The conditions obtained in the equivalence case are re-encountered 
now. Moreover, we obtain two additional, and very complicated, conditions on the 
functions gi{x) to have compatibility. One of these two conditions is the one satis- 
fied by system (IT^ as a particular case. We have seen that system ( IT^ cannot be 
seen as a particular case of a linearizable family of systems (fTSll by the equivalence 
method. 

Example 6. We consider a planar differential system of the form: 

i = 9oix) + giix)y + g2ix)y^, y = g3{x)y + g4{x)y'^ , (16) 

where gi{x), i = 0,1,2,3,4, are arbitrary functions. We remark that this system 
has y = as invariant algebraic curve and we propose an inverse integrating factor 
which contains this information. We are going to give conditions on the functions 
gi{x), i = 0, 1,2,3,4, such that the system has an inverse integrating factor of the 
form: 

V{x,y) = hi{x)y + h2{x)y^ + hs{x)y^ , 

where hi{x), z = 1, 2, 3 are suitable functions. The imposition for V to be an inverse 
integrating factor of system flTB]) gives rise to a polynomial in y of degree 5. From 
the coefficients of y^ and y^ of this polynomial we deduce that: 

hsix) = k2g2{x), hi{x) = hgoix), 

where ko and ^2 are arbitrary constants. We end up with only three conditions 
which involve the functions gi{x), i = 0,1, 2, 3, 4 and the constants ko,k2. Since the 
function h2{x) is not concerned, we take h2{x) = for simplicity. The vanishing of 
this three conditions gives the following planar differential system: 

X = k2gi{x){go{xy + go{x)gi{x)y + gi{x)^y^), ^^^^ 

y = yikogoix) - k2go{x) - k2ygi{x)){gQ{x)gi{x) - go{x)g[{x)), 

where fco, ^2 are real constants and go{x), gi{x) are analytic functions. This system 
has the inverse integrating factor V{x,y) = ygi{x){kogo{x)^ + k2gi{x)'^y'^). 

The following rational change of variables x ^ z with z = 1 — k2 + ygi{x)/go{x) 
transforms the orbital equation associated to the system to the following linear 
differential equation: 

dy ^ yz^ 

dz ~ (1 - k2 - z) {k2 + k^il - k2 - z)^y 
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2.4 Inverse integrating factors of only one variable 

Example 7. In the following example we describe another method of linearization, 
which consists on imposing conditions to the system so as to obtain an inverse 
integrating factor depending only on y. 
Let us consider the following system: 

X = y + Ci{x) + C2{x)y'^ , y = c^ix) + C4{x)y + C::,{x)y'^ , (18) 

where Ci{x), i = 1,2,3,4,5 are arbitrary functions. Let us impose a function 
V = Voly) to be an inverse integrating factor and we look for conditions on Cj(x), 
i = 1,2,3,4,5, to accomplish this fact. The condition for Vo{y) to be an inverse 
integrating factor for (fT8|) reads for: 

(c3(x) + C4{x)y + C5{x)y^) V^iy) = {c[{x) + c^ix) + 2 Cr,{x)y + ^(x)?/^) Vo{y). 

Imposing that Cj{x) = kj c[{x) for j = 3,4,5 and C2{x) = ki + k2Ci{x) with ki, 
i = 1,2,3,4,5, arbitrary constants, we get a linear differential equation for Vo{y). 
We change Ci(x) to c{x) for simplicity in notations and we have that the system: 

X = c{x) + y+{ki + k2c{x))y'^, y = c'(x)(A;3 + k^y + k^^y"^), (19) 

linearizes since it has an inverse integrating factor V = Vo{y) which needs to satisfy 
the following linear differential equation: (^3 + kiy + k^y'^)VQ{y) = (1 + A;4 + 2k^y + 
k2y'^)Vo{y). We note that the function Vo{y) is the exponential of the primitive of a 
rational function, that is, it is of Darboux type. 

We are going to present a change of variables to system (HM which transforms 
the corresponding orbital equation to a linear differential equation. The following 
rational change of variables x ^ z with z = c{x) + y + {ki + k2c{x))y'^ transforms 
the orbital equation associated to system (fT9|) into: 

dz _ 1 + 2kiy - k2y^ ^ 1 + 2/1:2^31/ + 2k2y\l + A;4 + k^y) + k^y^ ^ 
dy 1 + k2y'^ (1 + k2y'^) (k^ + k^y + k^y^) 

which is linear. Thus, we get that system (fT9l) is a particular case of the families 
described in [8]. 

It is evident that linearizable systems of the form (fT9!) of any degree in y, can 
be constructed in an analogous way. The rational change of variables which would 
transform it to an ordinary differential equation of linear type, and thus relate it to 
the work is x 1— >• 2; with z = P{x,y), where P{x,y) is the function defined by 
X = P{x, y). 

This example suggests the following question: does it always exists an algebraic 
change of variables that transforms a linearizable system with a Darboux inverse 
integrating factor of only one variable to a system whose orbital equation is linear? 



18 



2.5 Inverse integrating factors of the form V{x, y) = r{x)h{y). 



Example 8. The following examples of linearizable systems are obtained imposing 
that the inverse integrating factor is a product of two functions: one in the variable 
X and the other in the variable y. 
Let us consider the system 

X = gi{x)fi{y) + g2{x)f2{y), y = 93{x)f3{y) + 94{x)h{y), (20) 

where gi{x) and fi{y) are arbitrary functions, i = 1,2,3,4. We impose this system 
to have an inverse integrating factor of the form V{x,y) = r{x)h{y). When we 
substitute this expression in the partial differential equation which defines an inverse 
integrating factor, we get: 



My) + My) ^ ) h'{y) + 



f'^{y) ^ j^(^) f 92{x)r' {x) - g'^{x)r{xY . ^ 



gsix) V r{x)g3{x) 

gi{x)r'{x) — g'i{x)r{x) 



r{x)g3{x) 



h{y) = 0. 



+ h{y) 

We take 

\gi[x) / /Com \gi[x)/ 

where ki, i = 0,1,2 are real constants, and the previous relation reads for: 

My) + ^) h'{y) + (^koMy) - My) - ^) Hy) = o, 

which is a linear ordinary differential equation for h{y). We obtain the following 
system: 

X = kogi{x){gi{x)fi{y) + g2{x)f2{y)), ^^^^ 

y = ihMy) + My))igiix)g2{x) - gi{x)g2{x))/ki, 

where ko, ki are real numbers and gi{x) , g2{x) and fi{y), i = 1,2, 3, 4, are arbitrary 
functions. This system has the inverse integrating factor V{x,y) = gi{xYh{x), 
where h{x) satisfies the aforementioned linear differential equation of order 1. 

The following change of variables x ^ z with z = g2{x)/gi{x) transforms the 
orbital equation associated to system fl?Il) into the following linear differential equa- 
tion: 

^ ^ _^ My) + zMy) 
dy '^kiMy) + My)' 

As in the previous example we are addressed the question of the existence of an 
algebraic change of variables that transforms a linearizable system with a Darboux 
inverse integrating factor of the form V{x, y) = r{x)h{y) to a system whose orbital 
equation is linear. 
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3 More general inverse integrating factors 



In the first sections of tliis work we liave provided several examples to exhibit that 
certain families of systems ([1]), which are polynomial in the variable have an 
inverse integrating factor which is a polynomial in y. In this section we treat other 
expressions of an inverse integrating factor which contain the polynomials in y as 
a subclass. This generalization is done in two steps. We first introduce a real 
parameter a which allows to study polynomial inverse integrating factors which 
are a polynomial in |/ up to a real power. This real parameter does not involve 
any change in the linearization process described so forth. The second step is to 
consider inverse integrating factors which are a power series in y and we end up 
with a numerable set of linear differential-difference equations. 

We remark that the linearization process can also be carried out by imposing a 
function Va{x,y) of class in some open set U of M^, non locally null and which 
satisfies the following partial differential equation: 



dx dy \ dx dy 

where a is a real number. In the case a = 1 we recover the method of linearization 
described in the previous sections and this real free parameter a gives a general- 
ization of the linearization process which can lead to wider families of linearizable 
systems. The knowledge of a function Va{x,y) satisfying this partial differential 
equations gives that V = Va'^ is an inverse integrating factor of system ([T]). In the 
following, we do not impose a system to have an inverse integrating factor but to 
have a function V^(x, which satisfies the partial differential equation (!22|) where 
a is a real parameter. 

In order to illustrate this linearization process we describe an example, where 
the expression of the function V"a(x, ?/) is not a polynomial in y but a power series 
in that is, 

oo 
n=0 

Let us consider a quadratic polynomial differential system: 

X = aoo + cliqx + aoiy + 020^;^ + auxy + ao2l/^, ^^^^ 
y = 600 + hox + boiy + b2ox'^ + buxy + 6o2l/^ 

where a^j and bij are real numbers. We impose that Va{x,y) = J2n>o'^ri{x) y"', 
satisfies the corresponding partial differential equation fl22|) . We fix a natural num- 
ber n, we equate the coefficients of y"' in the development of (1221) and we get the 
following linear differential-difference equation for Vn{x). As in the rest of the work. 
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v'^{x) means the derivative of f„(x) with respect to x. 



(aoo + aiox + a2Qx'^) v'„{x) + (aoi + aux) + ao2 ^^^2(2^) 

+ {n + l) [boo + biox + h2QX^) Vn+i{x) + [(n - a) (601 + hix) 



(24) 



— Q;(aio + 2a2oa;)] + ((n — 2a — 1)602 — actii) "i^n-ila^) = 0. 

We include several examples of planar systems whose integrability can be deter- 
mined with a function V"q,(x, ?/) which is a power series in y. 

Example 9. We are going to take f„(x) in recurrence ^I^., of the form 
Vn{,x) = q{x) pn{x) where q{x) is a suitable function in x (usually of Darboux 
type) and Pn{x) is a polynomial in x. In particular, in this example we impose a 
function of the form f„(x) = q{x) ipn Hn{x) to be a solution of equation (l24l) . where 
q{x) is a suitable function, ipn is a suitable sequence of real numbers and Hn{x) 
is the Hermite orthogonal polynomial of degree n. For further information about 
orthogonal polynomials and the identities they satisfy, see for instance The 
Hermite orthogonal polynomials satisfy the following two identities: 



We substitute the expression f„(x) = q{x) ipn Hn{x) in and we use the previous 
identities to simplify it. We have that the following relation must be fulfilled: 

ao2q{x)ipn-2H'^-^2{^) + Ai{n, x) Hn-i{x) + ^2(72, x)if„_2(x) = 0, 

where the Ai{n,x) are expressions involving the sequence the function q{x) 
and the parameters aij, bij and a. Since we have already used all the identities 
of the Hermite polynomials, we need that the coefficients of H!^_2{x), Hn-i{x) and 
Hn-2{,x) independently vanish. We impose that ao2 = and the following values 
make that Aj(n, x) = for i = 1, 2: 



Hn+i{,x) = 2xHn{x) - 2nHn-i{x) 



H'^ix) = 2nH^^^{x). 



(25) 




601 



bii — —020, h 



0, a 



-1, 



'02 = u. 




with /i = 2 




+ 2 



(aii^oo ~ c^oi^io) 

•^20 



- 1. 
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We rename the free parameters by aoi = ^oan, 020 = — a2«ii, &00 = boau, 
bio = biaii, 620 = ^2Cfii and an is taken to be an = 1. We obtain that the 
quadratic system: 

X = (ao + x){y - a2x), y = Bq + bix + b2x'^ + 02(00 + x)y, (26) 

has the following expression Va{x,y) which satisfies the partial differential equation 
(122]) with a = -1: 

Va{x,y) = q{x) ^fnHn{x)y'^ = q{x)a2 ^ — Hn{x) ( — ) . 

^ — ^ ^ — ^ n\ V 02 / 

n>0 n>0 ^ ^ 

The Hermite polynomials have the following generating function: 

exp{2xy-i/2} = XI ^^«(^)^"' 

n>0 

which also appears in the book [Ij. This identity allows us to identify the power 
series given by V^(x, ?/) and we obtain the following inverse integrating factor V = 
for system fl26l) 

-a2— 2bQ+2QQb-[^— 2aQb2 

which is of Darboux type and which is not a polynomial in any of the variables x 
or y. 

In this example, we have solved recurrence fl24l) by imposing it to be compatible 
with the identities (12^ . Using the solution of the recurrence, we have encountered 
the inverse integrating factor of system fl2B|) and we have, therefore, integrated the 
system. 

Example 10. The system x = 1 — x^, y = |/(x — ?/), has an inverse integrating 
factor of the form V = Va'^ with a = — 1/2 and where: 

1 1 — xy 

K.(a;,i/) 



Vix^-)- exD { y'-^<'2xy-b2x' + 2{aob2-b,)x 



(1 -a;2)i/4 1 _ 2x7/ + 7/2- 

This function V^(x, has been determined using that it satisfies; 

00 

Va{,x,y) = q{x) Xr„(x)t/", 

n=0 

where q{x) = (1 — s^)"^/^ and Tn{x) is the Chebyshev polynomial of first kind and 
of degree n. This choice makes the recurrence equation ( l24l) to be satisfied in this 
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case. As stated in the book pQ , the Chebyshev polynomials of first kind satisfy the 
following identities: 

Tl 

Tn+i{x) = 2xTn{x) -Tn-i{x), r^(x) = 2nTn-i{x) + ^^^—^Tl^_2{x) , 



1 — xy 



,2 



Y.^n{x)y\ 



l-2xy + y 



Example 11. The system x = 1 — x^, y = 1 — xy, has an inverse integrating 
factor of the form: 

V[x,y) = , 
a/1 - 2xy + ?/2 

which has been determined using that it satisfies; 



V{x,y) = q{x) Y,Pn{x)y\ 



n=0 



where q{x) = (1 — x^)^ and Pn{x) is the Legendre polynomial of degree n. As stated 
in the book [1], the Legendre polynomials satisfy the following identities: 

{n + l)Pn+i{x) = {2n + l)xPn{x) - nPn^iix), 

Tl 

^«(^) = (xP„(x) - p„_i(x)) , 

^ OO 

Example 12. The system 

X = 1 + {2a^ - 7)?/^ + 6?/^ - 2axy{\ + y^) + 2x^y^, 
y = V(l 

where a is a real parameter has an inverse integrating factor of the form: 

1 ( 2axy - {a'^ + x^)y^ 

''(^•^' = 7r3^'='n 

This function V{x, y) has been determined using that it satisfies; 



(27) 



V{x,y) = y 



n\ 2-^ 

n=0 
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where Hn{x) is the Hermite polynomial of degree n. The identity 



1 




is called Mehler's-Hermite polynomial formula, see for instance [2]. 

These examples suggest that there is a connection between some Darboux inverse 
integrating factors and orthogonal polynomials. 
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